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Abstract 



We show how to generate coupled KdV hierarchies from Stackel separable systems of Benenti 
C ■ type. We further show that the solutions of these Stackel systems generate large class of finite-gap 

and rational solutions for cKdV hierarchies. Most of these solutions are new. 

■ 1 Introduction 

00 ■ 

CO ' In T we presented a systematic method of passing from Stackel separable systems to infinite hierarchies 

of commuting nonlinear evolutionary PDE's. We presented the idea in a concrete case of Stackel systems 
lO ' of Benenti type. In this paper we recognize the obtained hierarchies as the well known coupled Korteweg 

I - de Vries (cKdV) hierarchies of Antonowicz and Fordy [2] written in a different representation. We also 

clarify and significantly simplify the approach developed in [T] . The main idea of the paper is however to 
present a new method of generating solutions of soliton hierarchies from solutions of the related Stackel 
systems. 

From the very beginning of development of theory of integrable systems in the late 60 's major efforts 
}J] ' have been put into constructing various ways of finding their solutions. Among many others, a possible 

way of finding solutions of integrable systems is through various kinds of symmetry reduction, where 
one starts from an infinite-dimensional integrable system and obtains after such reduction an integrable 
ODE. If one then succeeds in solving this ODE (for example by finding separation coordinates in case of 
Hamiltonian systems) one can then reconstruct the corresponding particular solutions of the integrable 
PDE. This method originated in ^3 and has been developed in [5], [1], [5], [5], [7] (see also [8J) and later in 
a large number of papers. In this paper we present an opposite approach in the sense that we start with 
large classes of Stackel systems written in separation coordinates so that their solutions are explicitly 
known and in few steps we construct from these Stackel systems an infinite hierarchy of commuting 
evolutionary PDE's while the solutions of the considered Stackel systems become particular multi-time 
solutions of the systems of the obtained hierarchy. In this way we produce both well-known and new 
finite-gap type solutions of the KdV hierarchy as well as new rational and finite-gap solutions of cKdV 
hierarchies. In many cases this method also leads to implicit solutions of cKdV hierarchies. 

The paper is organized as follows. In Section 2 we briefly describe the starting point of our consid- 
erations, that is Stackel separable systems of Benenti type, including their general solution. In section 3 
we relate with our Stackel systems a class of weakly-nonlinear semi-Hamiltonian systems (i.e. a class of 
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hydrodynamic-type systems) that are reductions of the so caUed universal hierarchy [5] ■ These systems 
are in our formulation defined by Killing tensors of Stackel metrics. In section 4 we explicitly construct 
hierarchies of commuting evolutionary PDE's and present a transformation that maps these hierarchies 
onto the well known cKdV hierarchies of Antonowicz and Fordy [2] • This section contains the main result 
of this paper, i.e. Theorem[8l that produces large families of solutions to our coupled hierarchies. Finally, 
section 5 is devoted to studying specific classes of solutions that we call zero-energy solutions that contain 
both rational and implicit solutions of our hierarchies. 

2 Stackel systems of Benenti type 

Let us consider a set of canonical (Darboux) coordinates (A, /i) = {Xi . . . , , Xn, fJ-i, ■ ■ ■ , ^J■l) on a 2n- 
dimensional Poisson manifold M. A set of n relations of the form 

(/3i(Ai,^i,ai, . . . ,a„) = 0, i = 1, . . . ,n, ai e R (1) 

(each involving only one pair A^, fj,i of canonical coordinates) are called separation relations |10j provided 
that det ^ff^^ ^ 0. Resolving (locally) equations ([1]) with respect to we obtain 

ai^Hi{X,fi), i = l,...,n. (2) 

with some new functions (Hamiltonians) Hi{X, fi) that in turn define n canonical Hamiltonian systems 
on M : 

A, = ^, ^^U-~^, ^-h-,n (3) 

(here and in what follows the subscript denotes derivative with respect to the subscript variable). From 
this setting it follows immediately that the Hamiltonians Hi Poisson commute. The corresponding 
Hamilton- Jacobi equations for all Hamiltonians Hi are separable in the (A, /x)-variables since they are 
algebraically equivalent to the separation relations ^ . 

In this article we consider a special but important class of separation relations: 

n 

^a.Ar^' =^/(A.)m? + S7(A.), *-l,...,n, (4) 

where A and B are two constants. Notice that since the functions 7 and / do not depend on i the relations 
(|4]) can in fact be considered as n copies of a curve - the so called separation curve in A-^ plane. The 
Hamiltonian systems obtained from this class of separation relations have been widely studied and are 
also known as Benenti systems. Benenti systems constitute the simplest, but still very wide, class of all 
possible Stackel separable systems. It can be shown [11 that this class contains all quadratic in momenta 
Stackel separable systems since all other systems of this type are constructed from ^ by appropriate 
generalized Stackel transforms and related reciprocal transformations. Below we remind some established 
facts about Benenti systems. 

The relations ^ are linear in the coefficients a^. Solving these relations with respect to we obtain 

ai = A^i^K,G{f)ii + BVdn)^H,, z = l,...,n, (5) 

where we use the notation A = (Ai, . . . , A„)^ and fi — (/ii, . . . , /x„)^. Functions Hi defined as the right 
hand sides of the solution ([5]) can be (locally) interpreted as n quadratic in momenta /i Hamiltonians on 
the phase space M = T* Q cotangent to a Riemannian manifold Q equipped with the contravariant metric 
tensor G(/) depending on function / only. As mentioned above, these Hamiltonians are in involution 
with respect to the canonical Poisson bracket on T*Q. Moreover, they are separable in the sense of 
Hamilton- Jacobi theory since they by the very definition satisfy Stackel relations ((T|). The objects Ki in 
([5]) can be interpreted as (1, l)-type Killing tensors on Q related to the family of metrics G{f). The scalar 
functions Vi^"/) depend only on the function 7 and can be considered as separable potentials. Further, 
the metric tensor G and all the Killing tensors Ki are diagonal in A-variables. More specifically, in 
A-variables they attain the form 
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G(/)=diag(^,...,^) with A, = n (A. -A,) (6) 



Ai A„ 



and 

respectively. Here and below qi = qi{X) are Vietc polynomials (signed symmetric polynomials) in A: 

ft(A) = (-l)* ^ A,, ...A,., ^ = l,...,n. (7) 

l<si<S2<...<Si<n 

that can also be considered as new coordinates on the Riemannian manifold Q (wc will then refer to 
them as Viete coordinates). Notice that the Killing tensors do not depend on a particular choice of / 
and 7. 

Remark 1 The general n—time (simultaneous) solution for Hamilton equations associated with all 
the Hamiltonians (0i attains the form 



1 ^ \ f X^~^ 
c, = ±;;-7=y^ I I dXr, i = (8) 



To see this it is enough to integrate the related Hamilton- Jacobi problem. Now, with every Hamiltonian Hi 
in (l5|) we can associate the corresponding inverse Legendre mapping (fiber derivartive) : T*Q ^TQ 
given in the natural coordinates (A,/i) on T*Q and (A, Af.) on TQ respectively by 

Cr\\,fi) = {X,2AK,Gfi) - (A, A*,). 

Performing n corresponding Legendre transforms we obtain n Lagrangians Li : TQ — s-R given explicitly 

by 

Ui\,Xu) = -^XlgK-'Xu - BV,i-f), (9) 

where g = is the corresponding covariant metric tensor. These Lagrangians give rise to n systems 
of Euler-Lagrange equations 

Et,{L,) = 0, z = l,...,n (10) 

where Et^ is the Euler-Lagrange operator with respect to the independent variable ti. By construction, 
the solutions ([8]) are also general solutions for all the Euler-Lagrange equations pO|l . It means that for 
a particular i the general solution of Euler-Lagrange equation Et-{Li) = is given by ^ where tj is 
constant for j ^ i. 

3 Dispersionless Killing systems of Benenti type 

We have shown above that separation relations ([4]) lead to separable Liouville systems (Benenti systems) 
and we also presented its general solutions. It turns out that with the set of Killing tensors Ki of any 
Benenti system we can also relate a set of first order evolutionary PDE's. 

Definition 2 [12] For any fixed j G {1, n}, any of the following n systems of evolutionary PDE's of 
the form 

Xu^K,KJ^Xt^=Z,,{X,Xt^), 1 = 1,..., n. (11) 

(where X = (Ai, . . . Xn)^ ) with Ki being Killing tensors of a Benenti system ^D^, will he called a disper- 
sionless Killing system of Benenti type. 
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The chosen variable tj in pT|) plays the role of a space variable while the remaining variables ti 
should then be considered as evolution parameters (times). Equations pT|) constitute a set of n integrable 
dispersionless equations, belonging to the class of so-called weakly nonlinear hydrodynamic-type systems, 
i.e. that are semi-Hamiltonian in the sense of Tsarev [13], [14] and linearly degenerate [l^, where the 
variables Xi are the Riemann invariants for the system. For a chosen j, the systems (jlip can be considered 
as n vector fields Zij (with a fixed j) on some infinite dimensional function space Mj of functions 
(Ai(tj), . . . , Xn{tj)). It can be shown that 

Proposition 3 On Mj (i.e. for a fixed j) the vector fields commute: 

[Zij,Zkj]=0 i,k^l,...,n. 

The following proposition is crucial for our study: 

Proposition 4 Every mutual solution X(ti, . . . ,f„) (given by ^) of all Hamiltonian systems with 
Hamiltonians of Benenti type ^ is also a solution of all corresponding Killing system in ill]) . 

Proof. For the class of Benenti systems ([S]) the spatial part of ^ attains the form 

f) M 

Xu^^ = 2AK,G^i, i = l,...,n. (12) 

So, for any fixed j E {1, ...,n} we can eliminate the momenta /i from (H^. This yields pTjl . ■ 
Thus, all the solutions ([8]) are also solutions of all n Killing systems in (|lip . Moreover, we have 

Theorem 5 The n-time general solution of all the Killing systems in ill]) is given by 

^ /" * 

U+Ci = y^ — j—dXr, i = l,...,n (13) 

( where o,re arbitrary functions of one variable ) 

The proof of this statement can be found in [15]. Thus (fT3|) indeed encompasses all the solutions ([8]). 
We also see that any solution (fT3)) can be written in the form ([8]) with appropriately chosen /, 7, a^, 
A and B) so that on a given surface of fixed values of all (for example for zero-energy solutions, see 
below) with any such solution we can associate infinitely many corresponding Stackel systems ([3]) sharing 
the same solution. 

Consider now solutions ([8]) as a specific class of solutions of (fTT]l . Since this class of solutions - by 
construction - satisfies all the Euler-Lagrange equations (I10|l we can treat these equations as additional 
bonds that these solutions satisfies. We can therefore use these bonds to express some variables Xi by 
other As. Thus, within the class ([5]) of solutions of the Killing system (ITT|) we can perform a variable 
elimination (reparametrization) that turns (jlip into entirely new sets of evolutionary PDE's. Below we 
demonstrate that in some carefully chosen cases this reparametrization turns systems ()11[) into systems 
with dispersion (soliton hierarchies) with the solution ([8]) being also a solution of these new systems with 
dispersion. Specifically, we will produce this way all coupled KdV hierarchies as well as new interesting 
classes of their solutions: finite-gap and rational solutions. We will also demonstrate that our hierarchies 
indeed are the well known cKdV hierarchies obtained by Antonowicz and Fordy through the energy- 
dependent Schrodinger spectral problem [2]. 

From now on we will assume that 

/==A", 7 = A^ m,/ceZ (14) 

so that (m attains the form 

n 

J2ci,K~' =A>Hl^ + BX^, z = l,...,n. (15) 
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We will denote the metric associated with / = A™ through ^ by G^"^\ It can be shown that for 
m — 0, . . . ,n the metric G*^™-* is of zero curvature while the metric G*^""*"^^ has non-zero constant curvature. 
The separable potentials associated with 7 = A'"' will be denoted by V^'^K The family of separable 
potentials V^'^'' can be constructed recursively [iQ by 

V,^'+'^ ^ - qM'^ with V}'^ = ^„ (16) 

(k) . . . (k) 

where we put 1^ —0 for i < or i > n. The first potentials are trivial: = Si,n-k for k — 

0, 1. The first nontrivial potential is V"/"' = ~qi, For k > n the potentials v}''^ become 
complicated polynomial functions of A. 

The Lagrangians ([9]) for any specific choice of m and n are denoted as /^"''"■'^ so that L"'"''"'' = 

From now on we will choose the representation j = 1 in (jlip so that the variable ti plays the role of 
the space variable. We denote therefore this variable as x: ti ~ x. The case of higher j is not discussed 
in this paper. Thus, we consider the Killing systems of the form 

\t^=K^\^ = Zl'i\,X^), i = l,...,n. (17) 

where the new upper index n in the ith vector field denotes the number of its components; note also 
that the second lower index in Z is now always 1 and can therefore be omitted. Also, from now on we 
denote the Lagrangian 2^"'™''^ simply as L"'™'*^ (to avoid the unnecessary index) so that 

i"'™'' = ^A;^5^'")A, - \ * = l,...,n. (18) 

In order to perform the aforementioned elimination procedure we will first pass to Viete coordinates ([7]). 
The Killing systems ([TT]) are tensorial so in Viete coordinates they have the form qt^ = Kr{q)qx or, 
explicitly 

d 

-j^qj = i<lj+i-i)x + i<lj+i-k-i)x - qj+i-k-i {qk)x) = {Z'i [q]y i, j = 1, . . . ,n (19) 

* k=l 

where Qq, = as soon as a > n and (Z" [q\y denotes the j-ih component of the vector field Z" [q\ = 
Z^{q,qx) (here and below the symbol f[q] will denote a differential function of q, that is a function 
depending on q and a finite number of its derivatives). One can see that iZ'^[q]y = (Z;[q])' for all 

1, j = 1, . . . , n. Further, in Viete coordinates the Lagrangian psp takes the form 

i"'™'' = i"''"''(g,'Z.) = ^9j.9^™^fe-S^/'\ * = 1,...,«, (20) 

where gj™^ = y^(2«-m-»-j) ^^^^ Euler-La grange operator E't^in (flO)) will be simply denoted as E, so 
in Viete coordinates ^ 

E = {El, . . . En) , Ei = - — . 

oqi 



Theorem 6 Lagrangian i20\) satisfies the following symmetry relations: 
1 . for a = 1 , . . . , 71 — 1 

(L"''"''=) = (L"'™+"^^-") , i = a + l,...,n, (21) 

that can also be written as 

E, (L"^"'''-) E,+a (L"'™-"^'^+") , i^l, n - a. (22) 

2. 

=£;;+! (L"+l'"'2"+'^+2) , cr = l,...,n- 1, Z = cr + 1, n. (23) 
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The proof of this theorem can be found in [T] . This seemingly technical theorem guarantees that the 
form of Euler-Lagrange equations survives the passage from rt-component to (n + l)-component Killing 
system and hence it will be crucial for the construction of soliton hierarchies below. The index a will be 
related with the number of components of the obtained soliton systems. 

For our further considerations we will also need a hierarchy of infinite Killing systems 

d 

-rrij = iQ3+t-i)x + ^ {qk {qJ+^-k-l)^ - Qj+i-k-i {qk).^) = {Zf [9])^ j = 1, . . .00 (24) 

* k=\ 

that is known as the universal hierarchy and has been considered in [S]. 



4 Coupled KdV hierarchies and their solutions 

We now shortly remind the reader our specific elimination procedure from py that turns the dispersionless 
Killing systems (fT9l) into cKdV hierarchies. More specifically, we show how to produce s (with s € N) TV- 
component {N G N) commuting vector fields (evolutionary systems) by eliminating some variables from 
a set of Killing systems (|19l) with the help of Euler-Lagrange equations for an appropriate Lagrangian 
£n,m,k rpj^g crucial for this procedure is that if applied to s -f 1 instead of s it yields the same set of 
s commuting vector fields plus an extra vector field that commutes with the first s fields. That means 
that this procedure leads in fact to an infinite hierarchy of commuting vector fields in the sense that 
for arbitrary s we can produce first s vector fields from the same infinite sequence of commuting vector 
fields. Moreover it turns out that this way we produce vector fields with dispersion (soliton systems), 
namely well-known coupled KdV hierarchies of Antonowicz and Fordy (in a different parametrization) . 
Details are as follows. 

Firstly, we choose A — 1 and B = —1. This specific choice of A and B is introduced only for a smoother 
identification of our systems with the aforementioned cKdV hierarchies; the elimination procedure works 
otherwise for arbitrary values of A and B. Consider all N possible splittings 

N=^a + a with cr e {1, ... ,A^} and a e {0, . . . , TV - 1} . 

Every such splitting leads to a separate hierarchy. Consider also the Killing systems (fT9|) . written in a 
shorthand way as: 

qtr = [qi, ■■■,qn], r = 1, . . . , n (25) 

where q= (gi, . . .g„)^. 

Remark 7 The first s = n — N+1 equations in I125\) are such that their first N components coincide with 
the corresponding components of the infinite Killing hierarchy \2J^ . The remaining n — s equations in 
i25\) are incomplete with respect to the infinite hierarchy {2^^ since beginning with the flow s + 1 systems 
\24-^ contain at its first N components also the variables qn+i, • • • , ^n+Af-i- 

Let us now choose m — —a and fc = 2n -|- in (fT4|) so that / — A^" and 7 = \'2n+N consider 
the last n — N Euler-Lagrange equations associated with ijn.-a,2n+N ^ q^^^ ^^^^ show [IJ that they have 
the form: 

En+I (L— ".2"+^) = 2g„ + V'i-'AT+lbl, •••,9n-l] = 0, 

En+2 (L— "'^n+A^j = 2g„_i + g„_2] = 0, 

i?„ (L"-".2n+A^) ^ 2qN+l + ^i"^ [gi, qN] = 0. 

Due to their structure, equations ()26p can be explicitly solved with respect to the variables qN+i, ■ ■ ■ ,qn 
which yields qN+i, ■ ■ ■ ,qn as some differential functions of gi, . . . , g^r: 

qN+i = fi"^ [gi, • • ■ ,gAf] 

: (27) 
gn = fi"-N+i [qi,---,qN]- 
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Naturally, the solutions ^ (with our choice of / and 7) solve both (P5)l and (P^ . Thus, within the class 
(0) of solutions il3\). we can use the Euler-Lagrange equations (|26p or rather their solved form (|27p to 
successively express (eliminate) the variables qw+i, as differential functions of gi, . . . , gjv in (^5]) . 

Plugging l[27|) into l[25|) produces n vector fields with N ~ a + a components: 

_ ^ ^ -„,Ar,a ^^^^^^^^^^ ae{0,...,iV-l} (28) 

(with q — (gi, . . . jQn)'^) ■ The higher components of (|25p disappear after this elimination within our 
class ([S]) of solutions. Moreover, since the first s equations in ([^5]) are complete in the sense of Remark 
[7] it can be shown that 

Zr m^Z^ [q\ r = l,...,s, 



meaning that the first s — n~ N +1 equations in ((28|) do not depend on n. Observe also that equations 
do not depend on n. Actually, if n increases to n' the last n — N equations in (|26p with this new n' 
will remain unaltered. This means that we can repeat this elimination procedure by taking n' ^ n + I 
instead of n (so that s increases to s + 1 and k ~ 2n + N increases to 2{n + 1) + N — k + 2 while a and a 
are kept unaltered. This new procedure (with n' = n + I instead of n) will therefore lead to a sequence 
of s + 1 autonomous N — (a + a)-component systems in which the first s systems will coincide with the 
corresponding systems obtained from the previous procedure (with s). This way we can obtain arbitrary 
long sequences of the same infinite set of commuting vector fields (soliton hierarchy) : 

qt^=Z/ [q] r = l,2,...cx3, a e {0, . . . , iV - 1} . (29) 
The second index a in ([M)) denotes different hierarchies. It can be shown [1^ that vector fields Z^'"commute 



for any i, j E N. 



Note also that functions /■^"■' in (^7)) depend on a so that indeed the procedure leads to N different 
hierarchies. Now, the n functions Xi{ti, . . . , i„) given implicitly by the system of equations 

-. n ^ . n— 2+a/2 
1 ^ — V / Ar 

= 1 

with 



i, + c. = ±-2^/ 7=^dXr 1^1,. ..,n (30) 



are solutions of the first n equations of the A^-component hierarchy ([29]) with N — a + a. The reason 
is that equations ([50]) are just equations ([5]) with / = A~" and 7 = X^n+N they clearly solve all 
equations (j28p . Moreover, it can be shown, that these solutions are zero on <?„+i, . . . , qn+N-i expressed 
as differential functions of gi, . . . , g^r through an appropriate system (|26p (with n' = n + N —1). It means 
that ([30|) indeed solve the first n equations in ([29]) . 

Consider now the following infinite multi-Lagrangian "ladder" of Euler-Lagrange equations of the 
form 

El (i",™+J-i,fc-J+i) = E2 (L"^"+^'"2''''~^'+2) = • • • = ^^(^n,™+j-«-fe-j+«) (31) 

with fixed m,k E Zi and with j ^ . . . , —1, 0,1,... (the multi-Lagrangian form of (|3ip is due to Theorem 
HI). The equations ^ that we use for variable elimination are then a part of this infinite ladder with 
m = —a and k = 2n -\- a + a and with j = 1, . . . , n. The equations (|26p are the only equations in the 
ladder ([3T|l (for this specific choice of m and k) that allow for the elimination described above. All other 
equations are complicated polynomial differential equations with no obvious structure that do not allow 
for any elimination procedure. As a consequence, there exists more solutions of the type (j30p associated 
with all the Lagrangians for (3 = 1, ... ,n — 1. However, one can show that for (i = 2 the 

variable qn+N~i is not zero on these solutions so that these relations solve only first n — 1 equations in 
([29|) . More generally, for any (3 > 1 the obtained solutions Ai(ii, . . . , i„) will only satisfy first n — /3 + 1 
equations of the hierarchy ([29]) . We can thus formulate the following theorem. 



7 



Theorem 8 For any f3 £ {0, . . . , n — 1} and any N ~ a -\- a < n the n functions A^(ii, . . . , t^j) given 
implicitly by the system of equations 



C-i — 2 



/ d\r z = l,...,n (32) 

and with A^"'^''''' given by 

are solutions of the first n — f3 + I (all n for f3 — 0,1) equations of the N -component hierarchy 



The variables ti = a;, t2, ■ • ■ , in-/3+i in (|32p are "dynamical times" (evolution parameters) of the 
hierarchy ([29]) while the variables t„_^+2, ■ ■ ■ ,tn are just free parameters (i.e. the solutions (|32|1 do not 
solve flows higher than the flow number n — f3 + 1). Note also that, due to the structure of ([7]), all n 
functions Xi{ti, . . . ,t„) obtained in (f32|) are necessary in order to compute N functions qi{ti, . . . ,tn-p+i) 
that solve ([^^ . In the case = 1 the solutions ([5^ are finite-gap solutions for the KdV equation with 
the parameters Ei playing the role of endpoints of forbidden zones. For f3 > 1 these solutions are up to 
our knowledge new. For iV > 1 all the solutions ([32| are new. 



Remark 9 For a fixed a e {0, . . . , N — 1}, the following map 

Qy(N.2N) 



dqN+ 



l-r 



l,...,N-a (33) 



(where V^^^'^^-* denotes the separable potential Vi'^'^^ in the dimension N) transforms the hierarchy (2. 
to the hierarchy generated by the spectral problem 

(A"a2 + VI. = A^*. (34) 

This is the well known spectral problem of Antonowicz and Fordy leading to N -component cKdV hierar- 
chies, one for each a G {0, . . . , iV — 1}. Thus, N hierarchies I129\) are nothing else as (reparametrized) N 
cKdV hierarchies from J^. 



5 Zero-energy solutions 

Let us now investigate - on few chosen examples - the nature of solutions (j32p in the case that all Oi 
vanish (zero-energy solutions). In this case the solutions ([5^ can easily be integrated yielding 

1 " 

t,+c, = ±——. ^Ai-'-'^/^ (35) 

r=l 

and contain therefore no /3 (are the same for all /3 = 0, . . . , n — 1) and no N except in n s + — 1. 
Therefore, we obtain the following corollary: 

Corollary 10 For any N — a -\- a < n, the n functions Xi{ti, . . . ,t„) given implicitly by the system of 
equations 135\) are solutions of the first n equations of the N -component cKdV hierarchy 129\} . 

Naturally, the solutions ([55)1 also solve (on the surface Hi = for all i) all Stackel systems ([5]) for 
which /(A)7(A) = A^"+°'. Observe also that the solutions ([55)1 solve all the Euler-Lagrange equations for 
the infinite sequence of Lagrangians jj^^-o^+i^'^^+N-j where j £ Z. As a consequence, the second part of 
the map ((33)) is zero on solutions ((35)) . The reason for this is that due to Theorem [6] the expressions 

(L^'^-"-2^) , r^N-a + l,...,N 
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can be written as 

Er (L"^"'2n+a^ ^ r = n-a + l,...,n 
whereas all the above expressions are members of the sequence 7j"^-"+J.2n+w-j^ j ^ 2^ \\^qX: they are 
zero on solutions ([55)1 . This implies that as soon as a > the solutions ([55)1 in the representation of 
Antonowicz and Fordy, reduce to the solutions of the corresponding hierarchy with the same a but with 
a = 0. 

Example 1. We start with the case = 1. We wish to obtain the first s = 3 flows in ([29|l . There is 
now only one splitting N = a + a possible, namely cr = 1 and a = 0. This choice leads to the usual KdV 
hierarchy. We have to take n = s + iV— 1 = 3. The Killing systems ([25|) have in this case the form 



The Lagrangian L" -a.2n+<T 




- 9291,. 1 = Zl (36) 
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qiq3,x - q?,qi,x | zl 



l3,0J = 1 (-^2 _ ^2 ^ _ ig^g^^^(j2a: + ^91,2^93,2; + ^92,x + 29293 " ^qUa ^ ^qiql + 49? 92 ^ 9l 

and the Euler-Lagrange equations (|26p attain the form 

E2 (L3.o,7) ^ 2(?3 - 6gig2 + 4(75* + jqh + \qiqi.xx - \q2.xx = 0, 
E:, (£3.0,7) ^ 2(72 - Ml - kqi,xx = 0. 
Due to their structure, these equations can be solved with respect to 52, 93 yielding (P7)) of the form 

q2 — \qi,xx + ^qi , 93 = tk^l.^jxa;. + |9l9l,a:a: + |9i,k + |9l- (37) 

Substituting ([57)1 to the Killing systems yields the three one-component flows 

9i,ti — 9l,x — 1 

9i,t2 = -^qi,xxx + 39i9i,K = ^2' (38) 
1 5 5 15 2 -^1.0 

9l,t3 = -^qi,xxxxx + 29l,x9l,xa: + ^9l9l,a:a;a: + -^9l9l,a: = -^3 

which are just the first three flows of the KdV hierarchy. By taking larger s we can produce an arbitrary 
ntimber of flows from the KdV hierarchy. Now, according to Corollary [101 the formula ([35| yields some 
specific solutions of all three flows in ([55)1 . Explicitly, this formula reads (with x = ti, c,j = and + sign 
in ([23): 

X = -ELi^' = -Pi 

t2 = -^Eti^' - {pI ~ 3pip2 + 3p3) (39) 

^3 = -^ELl^i^ = -\ {Pl - 5(PlP2 - P3){pI - P2)) 

where zi = Aj i = 1,2,3 and where pi = Ei=i-^i' P2 = ^1^2 + ^1^3 + ^2^3 and p3 = ^12:223 are 
elementary symmetric polynomials in z^. The right hand sides of (j39p follow from Newton formulas: 

Er=i^r= E (-i)-+-+-+--m^^i±^^±:::±^^^prpr---Pn", m = i,...,n 

Q:l+2a2 4-.■■+^^Q^Tl— 

(40) 



9 



that can be easily extended to the case to > n by taking larger n and putting all higher pi equal to zero. 
The system can be solved explicitly yielding the 3-time solutions 



Pi = -X 

15i3 + 2x^ - 15a;2t2 

5(^3-3^2) 



P2 ^ rFZ^ — — (41) 



P3- 



15 {x^ - 3t2 
On the other hand, it is easy to see that 



,i=('A'-2f^. (42) 

VP3/ P3 

Plugging (dH) into ([1^ we finally obtain a 3-time solution of the first three flows of the KdV hierarchy. 
It has a rather complicated, rational form: 

, -3 (675^1 - 270^32;^ + 2a;" + Ql^xHl - 1350a;tl) , , 

910^.^2, ^3 = ^ 43 

' {-Ibhx^-Abtl + x^ + Abxhf 

By taking larger s we can in this way obtain s-time solutions of first s flows of the KdV hierarchy. Finally, 
the map p3p is in this case trivial and reads ui = 2qi. 

Example 2. Let us consider the two field case: N — 2. There are now two splittings possible: 
A^ = CT + a = 2 + and N = cr + a = 1 + 1. We consider only the first two flows s = 2 (i.e. only the first 
nontrivial flow) so that n — s + N— 1 = 3 as before and therefore the original Killing systems are as before, 
i.e ([36]) - we just consider the first two of them. Let us first consider the splitting A^ = CT + a = 2 + 0. 
The Lagrangian L".-",2n+<T+a ^ ^3.0,8 jg 



1 111 

+ ql- 6gi(72'?3 + 4^39? - + 6g|gi - 5g29i + 9i- 

(note that its kinetic energy part is the same as in L^.o,? a^]-,Qve). The formulas ([26]) contain only one 
equation that can be solved with respect to yielding (P7|l of the form 

3 1 

qs = 3(?ig2 - 2ql + -qi.xx- 

Substituting this to (|25p (with s — 2 now) yields the first two flows of the first (i.e with a = 0) 2- 
component cKdV hierarchy: 

d ( Qi \ _ ( qi,x \ _ 2^3,0 



''^ ^ 1^ J \ 1^-^ J ' , (44) 

_ 91 \ / q2,x \ ^^3,0 

dt2 \ q^ J y 2q2qi^x + 4:qiq2,x - Qqfqi.x + jqi,xxx J ^ 

The zero-energy solutions psp (again with all q = and with the plus sign only) attain now the form: 

V3 

X = -^2^i=iZi = -pi 

h = - jY^Li^i = {Pi - 3piP2 + 3p3) (45) 

^3 _ 1 /„5 r/„ „ „ \/„2 



QJ2i=iZi = (pI - HpiP2 ~ Pz){p\ - P2)) 
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where pi again denote elementary symmetric polynomials in Zj but where now zi = X^^. Solving (|45p 
yields the following 3-time solutions: 

pi — —2x 

P2 = 2x^ + 2t2 (46) 
4 c, 

P3 = --X-'-4:Xt2-2t3 

where the variable plays the role of a free parameter for eqs. ([44|) . Moreover, we have 

qi = , q2 = — . (47) 

P3 P3 

Plugging into ([Tf|) we obtain the following solutions for 

gi(x,t„^3)- 3^^^^^3 + 6^,^ ^ g2(x,^.,t3)^ 3^^^^^^3^g^^^ (48) 

Note that Corollary [10] implies that the functions (I48|) solve first n = 3 flows of the hierarchy ((29|) with 
N = 2 and a = 0. In order to compute this third flow we just need to take s = 3 in the elimination 
procedure. The result is 

,xxx 



dt3 \ q2 J V ^iii2qi,x - iSgi^i.j; + 6qtq2^x + jqiqi,xxx + 3q2q2,x + \q2, ' ^ 



Let us finally pass to Antonowicz-Fordy variables. The map (j33l) attains the form 

ui = 2qi, U2 = 2q2 - iqj 
and it transforms both systems in (|44| to the representation of Antonowicz and Fordy. Explicitly: 

(49) 



d I "1 \ _ I "l.a; 
dti 



U^ \ / U-i ^ \ —3,0 

^1 m 

U2 J \ U2.X • 



dt2 \ U2 J \ U2Ul,x + \uiU2^x + ' " 



In the u-variables the solutions (|48)) yield solutions for (|49)l and attain the form 

, , 6(^2+^2) / . ?>{Qxt3 ~ ^x'^ - QxH2 - ^tl 

ui{x,t2,h)^— — ^ , W2(a;,i2,t3' - 



3i3 + 2x3 + 6xt2 ' ' ^' (3t3 + 2x3 + 6xi2)' ■ 

Example 3. Let us now consider the case N ^ a + a ~ 1 + 1. Again, we look for the first s = 2 
flows in (|29|) with = 2, a = 1. We have to take n = s + Af — 1 = 3 and the rather lenghty Lagrangian 

Tn,-a.2n + a+a _ r3,-l,8 _ ^ ( ^fi On n ^ n \ n"^ -4- ^ ( n'^ rr \ rr n 

^ — - J (9l ~ ^9192 + 93) ^l^x + 2 l*?! ^ 92j 9l,2;g2,a; 

1 1 2 1 

- -^qiqi,xq3,x - ^<7i92,x + -^q2,xq3,x 

+ ql~ Gqiq2q3 + 4^39? ~ ql + 6(?29i ^ 5^291 + 9i 

(the potential part is of course the same as in L^^'-^'^ above). The elimination equations ([26| yield again 
only one equation that solved with respect to qs reads 

93 = ~ + 3gi'72 - jqiqi^xx + \q2,xx- 

Plugging this into two first flows in we obtain the first two flows of the second (i.e. with a — X) 
2-field cKdV hierarchy: 



i ^ , \ 31 

d I qi \ / q2.x \ ^ -^34 

''■^^ \ q2 J V 2g2'?i,a; + 491^2,2: - ^^i.^^i.^x ~ Qqlqi,x - lqiqi,xxx + jq2,r ' ^ 
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Now, the solutions ([55]) in this case attain precisely the form or (|¥T|) in solved form, since both n 
and (7 arc the same in both cases. However, (|5Qp are two-component, so in this case we need to express 
both qi and 92 as functions of pi : 

.i=f^)'-2^, 9.= f^)'-2^. (51) 
VPs/ P3 VPs/ P3 

Substituting (|4T|l into ((5T|l we finally arrive at a 3-time solution of (|50p with as a free parameter: 



q2(x,t2,t3 



-3 (675t| - 270<3a;^ + 2a;i° + 675a;4i| - 1350a;i^) 

4:5{x^ - 3t2){x^ + 15x% - 30^3) 
(-15*20:3 - 45ti + a;6 + 45^3)^ ' 



(52) 



As in the previous example, the above solutions solve first n ~ 3 flows of this cKdV hierarchy which 
means that they also solve the next flow in the hierarchy (with the dynamical time ^3). We will however 
not write it here. Finally, let us pass to the Antonowicz-Fordy representation. The map (|33p is now 



Ui = 2qi, U2 -- 

so that ([50)1 in Antonowicz-Fordy variables reads 



2^2 - 3ql - ^qi,xx 



dti 



d 
dtn 



Ui 
U2 



Ml \ / Ml,x 
"2 y ~ V "2,x 

U2,x + luiUi,^ -I- iUi 
1 



= z. 



^3,1 r 



(53) 



In the u-variables the solutions ([52| yield solutions for ([53]) and attain the form 

'6 (675^^ - 270^3x5 + 2x^° + QlhxHl - 1350a;t3) 



Ui{x,t2,h) 



Ahxh) 



U2{x,t2,t3) = 



which is nothing else than the solution of the first three flows of the KdV hierarchy (with = 1), 
in accordance with the observation at the beginning of this section. 

Remark 11 It can be shown that our method yields rational solutions only for (7 = 1 and a = 2. For 
a > 2 our method leads to new implicit solutions of our cKdV hierarchies . 

Example 4. Let us thus finally investigate the case A^ = 3 = it + q; = 3-|-0 that will lead to implicit 

^ have 

now the form 



(54) 



= 2 so that 


n = 


s + N- 


1 = 4. The first 




(q.\ 




( qi,x \ 




d 


92 




q2,x 


= Z\ 


dti 


93 




q3,x 






^ 94 ) 




\ qi,x ) 






(q.\ 




1 


q2.x \ 


d 


92 




q?.,x + 


91 92, 1 - 9291, 1 


dt2 


93 




qi,x + 


9i93,x - 9391, X 




^ 94 ) 




\ qiqi.x ~ q^q\,x ) 



z^ 



The Lagrangian L"' °'''^^+^ = ^^4,0,11 yigijg one elimination equation ([26|l : 



4,0,11 



) =2q4 + 12g29l ^ 69193 - 393 - 59i - -9l,a;:E 
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Solving this with respect to 54 we obtain 

94 = -6(7291 + 3^193 + ^(72 - + ^Qi^xx- 
Substituting this into ([54]) yields two first flows of the 3-coniponent cKdV hierarchy (|29p with a = 0. 

(55) 

93,x + 9192,2: - 929l,x I = 

29391.x + Miqs.x - Qqfq2.,x - 12919291,2; + 3^292,2 + I09?9i,x + \qi,xxx 

The solutions p5p are now (again with all Ci — and with the plus sign only): 



(56) 




X — 


2Z^i=l'^i 


4w- 


3/01P2 H 


^3p3) 




t2 = 






5(piP2 


-P3)(P?- 


P2)) 


<3 = 


y Z^i=l'^4 


-7 (^1 - 


7(PlP2 


-P3)((P? 


- P2)^ + P1P3) - 7p4(p? - 2pip2 + P3)) 


t4 = 


gZ^i=l'^4 


1 

Pg 

9 ^ 









■'^ - 


1 






"3 






1 






^5 


{pI 


.5 7 


1 






~7 


(pI 



— 1/2 

with Zi = Aj where Pg is a complicated polynomial of degree 9 in pi that can be obtained from Newton 
formulas This system can not be algebraically solved with respect to pi. However, if we embed the 

system (|56p in the system 

" -22t=i^i = ~Pi 

h = -gELi^j^ = -7 {Pi ~ 5(pip2 - P3)(Pi - P2)) + 5pip4 + 5p5 (57) 

+ 7p?P5 

where a is a parameter, and where Qg is a polynomial of degree 9 in pi, . . . , such that Qglp^^Q = Pg, 
then obviously the solution of (|57p with the condition ps = will yield the solution for ([55)1 . The system 
([57)1 is algebraically solvable and yields 

Pi ^ Ri{a,x,t2,tz,U) (58) 

where Pi are complicated rational functions of their arguments. The polynomial equation ps = yields 
then implicitly a (multivalued) function a = /(a;, ^2, ^3, ^4)- Now, using the fact that 

/P3V.„P2 /P2\ PlP3 ^ 2 (pf + 2p2 _ 4p2p2 + 4piP3 - 4p4) 

9i = - — + 2 — , 92 = — - 2^^ H , 93 - 



P4/ P4 VP4y P4 P4 Pi 



we arrive at an implicit solution of (|55p of the form 

P5{a,x,t2,h,ti) ^Q, qi^ri{a,x,t2,h,ti) , i = 1,2,3, 

that is thus determined up to the implicitly expressed function a — /(x, ^3, ^4)- The concrete formulas 
have been obtained with the help of Maple and are too complicated to present them here. The map 
does not simplify these solutions. 
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6 Conclusions 



In this article we presented a method of constructing coupled Korteweg - de Vries hierarchies from Benenti 
class of Stackel separable systems. Our method allows for producing certain class of solutions of these 
hierarchies from solutions of corresponding Stackel systems (Theorem |8] and Corollary [T0|) . For = 1 
and /3 = 0, 1 these solutions are known but for larger /3 and for larger N they all seem to be new. Also, 
for cr > 2 all our solutions are implicit in the sense described in Example 4 above. It has to be stressed 
that our method is general in the sense that other separation relations lead to other hierarchies like for 
example coupled Harry-Dym hierarchies. These possibilities will be studied in a separate paper. 

Considering the class of solutions constructed in the paper one arrives at natural question whether 
these solutins are somehow related with some class of symmetry reductions. Some hint is given in the case 
TV = 1 : the finite gap solutions ([32|) for /? = and /3 = 1 for the KdV can be obtained from its stationary 
flows constructed with the help of first two local Hamiltonian representations of the KdV hierarchy. This 
question is beyond the scope of this paper but it certainly deserves a separate study. 
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